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Artificially structured coatings are widely employed to minimize materials deterioration and corrosion, the 
annual direct cost of which is over 3% of the gross domestic product (GDP) for industrial countries. 
Manufacturing higher performance anticorrosive coatings is one of the most efficient approaches to reduce 
this loss. However, three-dimensional (3D) structure of coatings, which determines their performance, has 
not been investigated in detail. Here we present a quantitative nano-scale analysis of the 3D spatial structure 
of an anticorrosive aluminium epoxy barrier marine coating obtained by serial block-face scanning electron 
microscopy (SBFSEM) and ptychographic X-ray computed tomography (PXCT). We then use finite 
element simulations to demonstrate how percolation through this actual 3D structure impedes ion diffusion 
in the composite materials. We found the aluminium flakes align within 15° of the coating surface in the 
material, causing the perpendicular diffusion resistance of the coating to be substantially higher than the 
pure epoxy. 

Barrier coatings, an essential branch of the corrosion-resistant coating family, are deliberately structured to 
operate in extremely aggressive environments such as seawater 1 " 3 to prevent corrosion of substrate materi- 
als 4 ' 5 . Among these coatings, the aluminium flake pigmented epoxy coatings are widely applied to protect 
steel structures such as ocean-going vessels because of their excellent anticorrosive performance 1 , as determined 
by their diffusion barrier properties 2 . The shape and spatial distribution of the aluminium flakes, usually several 
microns wide and a few hundred nanometres thick, added to the materials significantly improve the barrier 
properties by decreasing the transport rate of corrosive substances such as ions, water (vapour) and oxygen 
through coating films 1,2,6 . Additionally, the aluminium can react with the hydroxide ions produced at the cathode 
in the corrosion cell, reducing the pH at the coating- steel interface and decreasing the rate of cathodic disbond- 
ment 7 . Meanwhile, the epoxy resin gives the coating film abrasion resistance and functions as a corrosion 
inhibitor for the aluminium flakes 8 . However, the 3D structure of the coatings 9,10 has not been accessed in detail 
before. 

Here we studied an anticorrosive aluminium epoxy barrier marine coating composed of aluminium flakes, talc 
fragments, iron oxide particles and epoxy resin (matrix material). The measured coating samples were prepared 
by, firstly, introducing the pigment ingredients including aluminium flakes, talc fragments and iron oxide 
particles into the epoxy resin using high shear mixing until the pigment particles are appropriately dispersed. 
The amine curing agent was then added into the pigmented epoxy resin at 1 : 1 stoichiometry. The epoxy resin and 
amine curing agent were mixed immediately prior to coating application. The coating material was then applied 
on the substrate and used to prepare two types of samples. The non-degassed coating was obtained by directly 
curing the applied material at 100°C for an hour and then curing at room temperature for one month. The 
degassed coating was treated by an extra degassing procedure in vacuum before experienced the same curing 
process as the non-degassed coating. The former procedure is the traditional method; the latter degassing 
procedure could lead to potential improvements of coating technology. 

Both aluminium flakes and talc fragments work as barriers when they align parallel to the coating surface. In 
commercial products, talc fragments substitute for aluminium flakes to reduce cost as well as the risk of electrical 
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sparking. The iron oxide particles are added to produce the desired 
colour and to improve solvent and chemical resistance of the pro- 
ducts 11 . Serial block- face scanning electron microscopy (SBFSEM) 12 
and ptychographic X-ray computed tomography (PXCT) 13 ' 14 were 
employed for the purpose, after we found that full-field transmission 
X-ray microscopy using both laboratory and synchrotron sources 
could not reveal the contrast between the various internal structures 
of the material. 

Results 

SBFSEM, a serial sectioning imaging method 121516 , was used to image 
a degassed sample. The technique was originally developed for 3D 
spatial structure investigation of biological specimens 17 " 19 , from 
ultrathin section sample preparation for transmission electron 
microscopy (TEM) 20 . It acquires a continuous series of parallel 
images of fresh surfaces created by mechanical cuts made by a dia- 
mond knife installed in an ultramicrotome 12 . Three-dimensional 
(3D) images of samples are then generated by stacking the corres- 
ponding series of images together. 

A degassed and a non-degassed sample of the same material were 
imaged by PXCT which combines scanning coherent X-ray diffract- 
ive imaging 21,22 with computed tomography 23 ' 24 . Both samples were 
first cut to size using a focused ion beam (FIB) tool. A tomographic 
projection at each rotation angle was reconstructed from multiple X- 
ray diffraction patterns using a difference-map phase-retrieval algo- 
rithm 22 . PXCT produces quantitative 3D images of customised field 
of view with a representation of the average electron density per 
voxel 13 ' 25 . It avoids mechanical sample destruction and possible cut- 
ting artefacts which, in certain cases, can be detrimental to nano- and 
micro-scale sciences. 

Figure la is a contrast-reversed back scattering electron (BSE) 
image of the degassed aluminium epoxy barrier marine coating from 



the SBFSEM measurement. Flake-shaped dark objects, grey objects 
and small black dots in it were identified as aluminium flakes, talc 
fragments and iron oxide particles, respectively, as segmented by 
commercial software "Avizo" and shown in figure lb as dark-grey 
objects, light-grey objects and bright small spots. The remaining 
parts of figures la and lb were attributed to epoxy resin. This seg- 
mentation was validated by the quantitative values of electron den- 
sity obtained with PXCT, as discussed below. The 3D spatial 
arrangement of aluminium flakes within the material is illustrated 
in figure le (and as the supplementary movie 1), in which we can see 
a parallel arrangement of the aluminium flakes that cannot be 
observed in 2D slice images. 

Figure 2a presents a single reconstructed tomogram slice of a 
different piece of the same degassed sample obtained by the PXCT 
measurement. Here, the electron density is determined accurately 
enough for the components to be easily identified by the respective 
peaks of the density histogram, shown in figure 2b. The long, slim 
grey objects with electron density 0.793 X 10 30 /m 3 are identified as 
aluminium flakes, the slightly brighter grey objects with electron 
density 0.905 X 10 30 /m 3 are talc fragments, the small dots are iron 
oxide particles, the other grey parts are epoxy resin and the black 
regions are surrounding air. Upon conversion 25 , the aluminium and 
talc mass densities are measured to be 2.73 X 10 3 kg/m 3 and 3.00 X 
10 3 kg/m 3 , compared with 2.70 X 10 3 kg/m 3 and 2.78 X 10 3 kg/m 3 
from standard tables. The high contrast and boundary well-defined 
features generated in the tomogram slices enable their direct 3D 
segmentation into spatial volume regions. Figures 2d and 2e are 
3D renderings (and as the supplementary movie 2) after segmenta- 
tion, analogous to figures Id and le, again presenting the parallel 
alignment of the aluminium flakes distinctly. PXCT demonstrates 
here its combined advantages of ptychography 21 , a phase- sensitive 
imaging technique, and X-ray computed tomography. 




Aluminium I I Talc I I Iron Oxide 

Figure 1 | Results for the degassed aluminium epoxy marine coating imaged by SBFSEM. (a) An SBFSEM single slice, (b) Segmented result from 
figure la. (c & d) Rendering of 3D spatial structure image of the coating, (e) A 3D image of just the aluminium flakes in the coating, (f) Quantitative 
orientation analysis of the aluminium flakes in the coating material. Deviation projections of individual flakes are plotted stereographically around the 
surface direction (centre). Dash circles denote net 15°, 30°, 45° and 60° deviations from parallel alignment to the surface. The scale bars in a and b are 
5 microns. 
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Figure 2 | Results for the degassed coating imaged by PXCT. (a) A single tomogram slice perpendicular to the rotation axis. Arrows point to some of 
the identified aluminium flakes in red and talc fragments in green. The scale bar is 5 microns, (b) Histogram of electron density distribution of the 
tomogram slice in figure 2a. (c) Quantitative orientation analysis of the aluminium flakes, plotted in the same way as figure If. (d) Rendering of 3D spatial 
structure image of the sample, (e) 3D spatial arrangement image of just the aluminium flakes in the material. 



According to the sample mounting geometry (see Methods), the 
xy plane in figures lc, d and e is parallel to the sample surface. The 
coating-steel interface in figure 2a is labelled and is assumed to be 
parallel to the sample surface. The surface normal direction of each 
aluminium flake is defined as the orientation of its major principal 
inertia axis which is along the shortest geometrical axis of the object. 
We find that the aluminium flakes arrange nearly parallel to the 
coating sample surfaces, mostly within 15° deviation, as shown in 
the orientation stereograms in figures If and 2c (see details in sup- 
plementary methods). A group of misoriented flakes, having around 
40° deviation, is highlighted by a green circle in figure 2c. These flakes 
form a physical cluster lying next to each other in the material, as can 
be observed in the top-right of figure 2a or the lower-right part of 
figure 2e (circled by dashed yellow ellipses). 

In addition, the lengths and volumes of individual objects in the 
materials can be obtained from morphometric analysis after the 3D 
image segmentation. Thus, the size of each object and the pigment 
volume concentration (PVC), one of the most crucial formulating 
parameters of coating products, of the applied and in-service coat- 
ings can be revealed. These parameters of the aluminium flakes are 
tabulated in the supplementary materials. 

An interesting discovery in the non-degassed aluminium epoxy 
marine coating is that the aluminium flakes tend to follow the edges 
of nearby air pores, as seen in figure 3. This correlation in structure 
has not been observed before, and it is different from the common 
expectation that the aluminium flakes should be parallel to the 
external coating surface. This may be due to a surface alignment 
effect at the edges of air pores which is similar to that at the external 
coating surface, which is caused by solvent evaporation during the 
drying process 11 . The volume exclusion effect may arise simply 
because when, the air bubbles in the coating material nucleate and 
grow, they push the aluminium flakes to their periphery. This co- 
localised structure may be responsible for the strong barrier function 
against corrosive substances in the pores, which are presumably the 
penetration point for external agents and which provide space for 
accumulation of corrosive substances. However, this structure has 



the potential hazard that, once the aluminium flakes become eroded, 
the anticorrosive performance/protective capability of the coating 
would suddenly diminish because of decomposition of the barriers. 

Simulations of the percolation properties of the coating were per- 
formed using COMSOL "MultiPhysics" software package, based on 
the actual, segmented 3D spatial structures instead of modelling 
structures 26 . We determined the ion diffusion process in the coating 
films, which is regarded as the rate- determining process governing 
the cathodic delamination of coatings from steel and its corrosion 2 . 
The simulations were done by calculating electric conductivity a in 
the materials following Ohm's law 'J = — a VV which has the same 
mathematical form as Fick's first law 'J = — D VC that can describe 
diffusion in 3D systems when the concentration of the diffusing 
species can be assumed to remain constant. J is the current density 
or (ion) diffusion flux, VV is the potential gradient (locally the 




Figure 3 | Results for the non-degassed aluminium epoxy marine coating 
from PXCT measurement, (a & b) Tomogram slices perpendicular and 
parallel to the rotation axis through the reconstructed 3D volume, 
respectively. In figure 3b, the image plane is parallel to the coating surface, 
(c) Rendering of 3D spatial arrangement of aluminium flakes and air pores 
in the sample. The scale bars in a and b are 5 microns. 
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Figure 4 | Finite element simulation results based on the actual 3D spatial structures of the coating measured in this study, (a) Simulated flow along the 
direction perpendicular to the coating surface presented as streamlines within the real structure obtained from the SBFSEM measurement, (b) A slice cut 
of the current density map generated from simulation in panel a. (c & d) Simulated flow along the directions parallel and perpendicular to the coating 
surface of the structure obtained by PXCT. (e) Rendering of the 3D volume of the coating structure used for simulations in panels (c) and (d). 



potential difference, V, divided by the conductor length), D is the 
diffusion coefficient and VC is the concentration gradient. The ion 
diffusion and electric conduction are also relevant because both ions 
and electrons are charged particles and the potential difference is a 
driving force for ion transport during the electrochemical process of 
steel corrosion. The electric conductivities were calculated along the 
directions perpendicular and parallel to the coating surface direc- 
tions to determine the material's diffusion coefficient anisotropy 

R ^perpendicular/ ^parallel* 

Figures 4a, c and d present the simulation results as streamlines of 
flow, which are the paths followed by ions from one side to the other 
of the coating films during diffusion. The ions flow along surfaces 
of objects and go around the obstructions in the structures, which 
clearly reveals the barrier properties of aluminium flakes against 
ions. This validates, with direct evidence, the pictorial matrix models 
of barrier coatings found in the literature 27 . The current density map 
(figure 4b) illustrates the ions have much higher concentrations at 
the edges of flakes which could potentially become the points that 
corrosion starts. It may be worth applying extra protective treatment 
to the edges of the aluminium flakes. 

The simulation based on the SBFSEM structure (figure 4a), 
6 microns thick, shows the diffusion coefficient anisotropy to be 
R = 0.516. The 10 microns thick structure obtained from PXCT 
(figure 4e) gave R = 0.358. The smaller anisotropy ratio acquired 
in the latter case could be caused by a single big aluminium flake in 
the structure (see figure 4e) which occupies the majority area of the 
cross-section and blocks the transport routes of ions more drastic- 
ally. Therefore, trying to avoid breaking aluminium flakes either 
during manufacture or application process would contribute to 
improving the final performance of the coating films. Considering 
the not-quite parallel alignment of the aluminium flakes to the coat- 
ing surface, we could expect that pure epoxy would give higher ion 
diffusion coefficient and smaller anisotropy ratio value against the 
surface perpendicular direction. This illustrates that the diffusion 
rate of ions across the aluminium flakes pigmented coating is heavily 
reduced compared with pure epoxy. 

Discussion 

The morphological analysis of 3D spatial structures and the simula- 
tions performed on them reveal the relations between the structure of 
coatings and their properties. The results from SBFSEM and PXCT 
validated each other and both methods were proved to be effective 
ways to measure 3D spatial structures of complex industrial 



specimens. Compared with PXCT, SBFSEM provides higher resolu- 
tion and avoids employing iterative phase-retrieval algorithms. 
However, it is a destructive method because the diamond knife has 
direct physical contact with the imaging surfaces, which could poten- 
tially lead to deformation of shapes and movement of positions of 
features in the samples. PXCT is found to be especially powerful at 
quantitative (structural) characterisation as it shows the ability to 
distinguish the four components in the material quantitatively by 
electron density. In principle, it is a non- destructive method which 
should preserve the sample structure well during measurement. By 
taking measured 3D structures into detailed simulations, we have 
demonstrated an approach for validating mechanistic assumptions 
and potentially provide a practical method to engineer the efficacy of 
anti-corrosion coatings by modelling electrochemical process in 
materials based on the actual structures. This could be used to evalu- 
ate the performance of materials, and shorten the product research 
and development lifetime, compared with carrying out endurance 
studies of original specimens, which can take years. We found that 
the degassed specimens have more compact structure comparing 
with the non-degassed one prepared by conventional coating tech- 
nology, which means the former would provide better anti-corrosive 
properties. The degassing treatment could be used to enhance coat- 
ing performance and thus to further reduce the risk of steel construc- 
tion failures such as oil tanker break-ups. 

Methods 

SBFSEM measurement. The degassed aluminium epoxy barrier marine coating film 
was mounted on an aluminium specimen pin (Gatan, Inc.) using cyanoacrylate glue. 
The coating surface was placed perpendicular to the axis of the pin. The sample was 
then trimmed down by a microtome to create a block face sizes approximately 
500 urn X 500 urn. After that, it was measured by the SBFSEM system, installed in an 
environmental SEM and operated at 2.5 kV with 0.4 Torr chamber pressure in a water 
vapour environment. The BSE signals from the sample were acquired with 10 |is 
dwell time per pixel of size 30 nm X 30 nm, at a working distance of 5.8 mm. The 
sectioning thickness was 30 nm per slice with a cutting speed 22 s per slice. 
Altogether 200 slice images (6 urn thick) with a field of view 30.7 um X 30.7 um were 
obtained. 

PXCT measurement. PXCT experiments were performed at the cSAXS beamline 
(X12SA) of the Swiss Light Source, Paul Scherrer Institut, using X-rays of photon 
energy 6.2 keV. FIB fabricated cubic-shaped coating samples were measured. X-ray 
diffraction patterns in the transmission geometry were recorded by a PILATUS 28 
photon-counting detector (pixel size 172 um X 172 um), placed 7.2 meters 
downstream of the samples. During the measurements, the samples were rotated 180° 
with 0.5° intervals, obtaining projections of the samples at 360 orientations. At each 
orientation angle, the samples were scanned along circular shells with a 1.4 um radial 
step 29 in a plane perpendicular to the X-ray beam to cover a field of view extending 
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beyond the sample edges. The projections were reconstructed by the difference map 
algorithm 22 , then post-processed, aligned, combined and reconstructed to form a 3D 
volume according to the procedure described in reference 14. 
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